Josephson current through a Kondo molecule 



Rosa Lopez/ Manh-Soo Choi,^'^ and Ramon Aguado'' 



^ Departament de Fisica, Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain 



o 
o 

(N 



(N 



o 



> 

in 

cn 
o 

o 



I 

o 
o 



J3 



'''Department of Physics, Korea University, Seoul 136-701, Korea 
^Department of Physics and Astronomy, University of Basel, KUngelbergstrasse 82, 4056 Basel, Switzerland 
''Teoria de la Materia Condensada, Instituto de Ciencia de Materiales de Madrid (CSIC) Cantohlanco, 2804-9 Madrid, Spain. 

(Dated: February 6, 2008) 

We investigate transport of Cooper pairs through a double quantum dot (DQD) in the Kondo 
regime and coupled to superconducting leads. Within the non-perturbative slave boson mean-field 
theory we evaluate the Josephson current for two different configurations, the DQD coupled m 
parallel and in series to the leads. We find striking differences between these configurations in the 
supercurrent as a function of the ratio t/T, where t is the interdot coupling and F is the coupling 
to the leads: the critical current 7c decreases monotonously with t/T for the parallel configuration 
whereas Jc exhibits a maximum at t/F = 1 in the serial case. These results demonstrate that a 
variation of the ratio t/T enables to control the flow of supercurrent through the Kondo resonance 
of the DQD. 

PACS numbers; 73.23.-b,72.15.Qm,74.45.+c 



Introduction. — A localized spin (magnetic impurity) 
forms a many-body singlet with the itinerant electron 
spins surrounding it, and hence is screened. This, known 
as Kondo effect, happens at temperatures T lower than 
the Kondo temperature Tr-, which is the binding energy 
of the many-body singlet state U • It manifests itself as a 
quasi-particle resonance, namely, a peak of width Tk at 
the local density of states (LDOS) at the impurity site. 
In mesoscopic structures, such as a quantum dot (QD) 
coupled to normal electrodes, it leads to the unitary con- 
ductance as predicted theoretically 2] and demonstrated 
experimentally Ever since, nanoscale systems have 
been of particular interest as they allow ultimate tun- 
ability and enable thorough investigations of the Kondo 
physics. 

It is more intriguing when one couples localized spins 
to superconducting reservoirs and brings about the 
competition between superconductivity and Kondo ef- 
fect 0. The subgap transport through the Andreev 
bound states 'of (ABSs) in nanostructures between 
superconducting contacts is strongly affected by Kondo 
physics, leading to the phase transition in the sign of 
Josephson current 0. The large Coulomb interaction 
prevents the tunneling of Cooper pairs into QD; electrons 
in each pair tunnel one by one via virtual processes 
Due to Fermi statistics, it results in a negative super- 
current (i.e., a vr-junction). However, this argument is 
only valid in the weak coupling limit, when the gap A 
is larger than Tk. In the opposite strong coupling limit 



{A/Tk ^ 1) the Kondo resonance restores the positive 
Josephson current Q- Lately, a S-QD-S system (S de- 
notes superconductor) has been experimentally realized, 
confirming some of these physical aspects 4] . 

The ongoing progress in semiconductor nanotechnol- 
ogy has brought more complicated nanostructures such 
as molecule-like double quantum dots 9] (DQD's). The 
latter is an unique tunable system for studying inter- 
actions between localized spins 10] . On one hand, 
they are basic building blocks for quantum information 
processing ll] as a source of entangled states On 
the other hand, they are a prototype of the tunable two- 
impurity Kondo model , as demonstrated experimen- 
tally 13,|lJ|. A natural step forward is to study the artifi- 
cial molecule coupled to superconductors. Choi, Bruder, 
and Loss studied the spin-dependent Josephson cur- 
rent through DQD in parallel in the Coulomb blockade 
(CB) regime and generation of entanglement between the 
dot spins. Here, we examine the strong coupling regime, 
both in parallel and in series. We find striking differ- 
ences between these configurations in the supercurrent 
as a function of the ratio t/T: the critical current Ic 
decreases monotonously with t/T for the parallel config- 
uration whereas Ic exhibits a maximum at t/T = 1 in the 
serial case. 

Model — A superconducting Kondo molecule is mod- 
elled as a two-impurity Anderson model where the nor- 
mal leads are replaced by standard BCS s-wave super- 
conductors. The Hamiltonian reads 
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Cfcctcr (c^icr) describes the electrons with momentum k and 
spin a in the lead a (the dot i); Uia- = d\^dia. (j>a i^a) 
denotes the superconducting phase (gap) for the lead a. 
t (Vika) is the tunneling amplitude between the dots (the 
dots and leads). For a DQD in series Vi,kR = V2,kL — 
with Vi^kL = V2,kR = Vq; for a DQD in parallel Vi^kR = 
Vi.kL = V2ML = V2.kR = Vq. are the single-particle 
energies on the dots tuned by gate voltages. Ui is the 
on-site Coulomb interaction on the ith dot, and U12 is 
the interdot Coulomb interaction. 

We will be interested in the limit Ui,U2 00, U12 ~ 0, 
and —Si ^ r so that (riicr) = 1 for each i — 1,2. 
The model in this limit is well described in the slave- 
boson language [l^ . We first write the physical fermionic 
operator as di^a = bifida, where fi,a{bl) is the pseud- 
ofermion(boson) which destroys(creates) one "occupied 
(empty) state" on the dot i. We then introduce two 
constraints which prevent double occupancy in each dot 
(Ui, U2 00) by means of Lagrange multipliers, Ai and 
A2. The resulting model is solved within the mean- field 
(MF) approach, namely, replacing bi{t) {bi) = hi and 
thereby neglecting charge fluctuations in the dots. This 
approach has been applied successfully to the Kondo 
molecules coupled to normal leads 10]. 

The MF version of J'lf is now quadratic and con- 
tains four parameters, i.e., 61^2 that renormalizes the 
tunneling amplitudes, Vi_ka = biVi^ka and t = tbib2 
and Ai(2) that renormalizes the energy levels ei,2cr = 
£i + Ai_2- They are determined from the solution of the 
MF equations in a self-consistently fashion. The MF 
equations become simpler in the Nambu-Keldysh space 

where = {cl^^,C-kc.a) and = (^fl,ha' 

the spinors for the conduction and localized electrons. 
From the equation-of-motion (EOM) of the boson fields 
and the constrains the MF equations read: 
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The system of MF equations l(2Jl is now written in terms 
of the 2x2 matrix lesser Green function (OF) for the dots 
Gi,ja{t-,t') = i{^]a{'t')-,'^i<j{'t)) and the lead-dot matrix 
lesser GF G,,fc„^(i,t') - i{'^l^„{t').^ra{t)) (with com- 
ponents = G, "electron-like" GF G^^ = is the 
anomalous propagator, G^^ = F\ and G^^ = G cor- 
responds to the "hole-like" propagator). Following the 
standard procedure, the lesser GF's are obtained apply- 
ing rules of analitical continuation along a complex time 
contour to the EOM of the time-ordered GF (for details, 
see Ref. [l^)- The two off-diagonal GF's, namely, Gi^j„ 
(with i 7^ j) lead-dot GF Gi^kaa can be cast in terms of 
the diagonal dot GF Gi^ia = Gia using the EOM tech- 
nique. By doing this Eq. Q becomes (for simplicity we 
omit the spin indices) 
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G'f(2)(e) + G'f(2)(e) 

^(ei(2) - £1(2)) + ^[Gf(2)(e)(e - ^1(2)) 

+G<2)(e)(e + ei(2))]=0, (3b) 

where = Eq + A^ are the renormalized levels and 
Ti — ifV the renormalized hybridization, which are equal 
to the Kondo temperature Tk for the coupled system 13] 
[r — 2ttpmV^ with pn being the normal-state DOS]. 
At equilibrium we can employ Gf — 2z/(e)ImG[ with 
/(e) being the Fermi function. G[ is determined by 
the lead-dot and the interdot tunneling self- 
energies. Thus, the matrix elements of the lead-dot self- 
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-zr>s(e)Aexp(z0„)/lelwithps(e) =je\/Ve^ -A^e{e~ 
A) as the superconducting DOS and = 2TT\Vika\'^PN- 
The interdot tunneling self-energy is ^^^^^2) = ^'^^z [32(1) ^ 

^2{i)a\^^' where ct^ is the z component of the Pauli 
matrices and gl is the matrix GF for an isolated QD. 
Since we deal with a symmetric structure Ai = A ^ = 
£0 and r = r^. Thus, for example for the serial configu- 
ration the diagonal components of the GF read 
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Hereforth, for simplicity we assume a symmetric struc- where 2?(e) = [A'^ — — Y^ — (eo ~ i)'^] — 4p(y^ 
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s^) with X = s(e) cos((?!)), Y = s(e) sin((?!)), A 
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FIG. 1: (Color online). Andreev bound states in units of 
A. Left panel: Parallel case for t/T = 0,0.25, 0.5, 1, 1.25, 1.5 
(from top to bottom). Right panel: Serial case for t/V = 
0.1, 0.25, 0.5, 1, 1.25, 1.5 (from top to bottom). In both cases, 
A = O.IT^. 



e [1 + s(e)/A], and s(e) = Ar/x/A2 - e2. The GF ^ de- 
scribes the discrete Andreev bound states in the subgap 
region (|e| < A) as well as the continuum spectrum above 
the gap (|e| > A). The Andreev states appear as poles of 
the GF; i.e., the solutions of I?(e) = (see Fig. QJ. Ac- 
cordingly, the Josephson current has two contributions, 
-^tot = ^dis + -^con; -^dis from the discrete Andreev states 
and /con from the continuum (see Figs. [3 and Ol. The two 
parts of the Josephson current for the DQD in series (see 
below and Ref. |l9j for the DQD in parallel) are given by 
(hereafter currents are expressed in units of 2|e|A/fi,): 



Icon = -2i2 sin(20)Im / de 



AV{e) 
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(5b) 



In Eq. (|5a|l the summation is over all Andreev states 
Ep S [— A,e_F] {ep is the Fermi energy). Interestingly, we 
will see below [Eqs. ^ and Q ] that in the deep Kondo 
limit (A <^ Tk), we recover the short-junction limit for 
the Josephson current through non-interacting resonance 
levels 20]. In this limit the continuum contribution is 
almost negligible. For A < Tk the contributions from the 
continuum part becomes considerable and the behavior 
deviates from the non-interacting case. 

Next we present our results for the DQD in parallel 
and in series, respectively. We will choose Eq — —3.25, 
D = 100 (bandwidth), ep = and different values 
for the rest of parameters. For these values = 
D exp — 7r|eo|/r = 0.0036. (All energies are given in units 
of F). 

DQD in parallel — The problem is greatly simpli- 
fied by the transformation Ckaa = {ck^a ± Ckacr)V^ and 
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FIG. 2: (Color online). t/V < 1. Left(Right) panel corre- 
sponds to the parallel(serial) case. [a(d)] discrete /dis, [b(e)] 
continuum /con and [c(f)] total supercurrent /tot versus (f) 
for t/T = 0,0.2,0.4,0.6,0.8,1 (t/V = 0.2,0.4,0.6,0.8,1) and 
A = O.lTj^. Currents are in units of 2|e|/A/i. The curves are 
arranged with increasing t/V in the direction of the arrows. 
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FIG. 3: (Color online). The same as Fig. HI except that t/V > 
1. 



/i(2)(T = (/ecr i/oCT)\/2- The MF Hamiltonian is mapped 
into two independent Josephson junctions ("even" and 
"odd" ) through effective resonant levels at Eq ± i. Each 
of the two resonant levels accommodates an Andreev 
state /i'e/o and carries Jose phs on current /tot(0) = 
— 2eA/?i(/e + lo) (see Ref [l^ for their explicit expres- 
sions). Typical profiles of Andreev states on DQD in 
parallel are shown in the left panels of Fig.Q] Josephson 
currents are shown in the left panels of Fig. 01 < F) 
and Fig. (i > F). We note that in the deep Kondo 
limit A ^ Tk the DQD in parallel behaves in effect 
like two non- interacting resonance levels 20]: The An- 
dreev states and the Josephson currents are given by 
the corresponding expressions for the short-junction limit 
= A [l-T,/<,sin2 (,/,)] 1/2 and 



Ie/o[<t>) = %/o sin 20 [l-%io siv? i 
4fV(£2 



]-l/2 



(6) 



respectively, with 7^/^ = 'i^ /{e'^/^ + F). For very small 
t/T both dots have their own Kondo resonances at ep 
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and the Josephson current resembles that of a ballistic 
junction. As t/T increases the even and odd Kondo reso- 
nances eo±? move away from and as a result Jtot {4>) di- 
minishes becoming more sinusoidal Ie/oi4') ~ ^/oSin2(/). 

Experimentally more accessible is the critical current 
Ic lU, which is given by /c(</.) » 2eA/ftS^g^^^„j[l-(l- 
7g)^/^] and hence decreases with t/T. Notice that in spite 
of the simple formal expression for I^/oy Ic depends on 
the many-body parameters, Tk and eo±t, in a nontrivial 
manner through the solution of Eq. ^ . 

DQD in series — A completely different physical sce- 
nario is found for the serial configuration. Here the 
even/odd channels arc no longer decoupled and cause 
novel interference. The manifestation of the interference 
can be first seen in the profiles of the Andreev states as 
depicted in Fig. (^). An important difference from the 
parallel case is the almost flat spectrum with values close 
to -E'e/o(0) ^ ^ (reflecting a very small supercurrent as 
seen below). As t/T increases the spectrum possesses 
larger amplitude, and for t/T 1 we eventually recover 
the spectrum of a ballistic junction Ee/o{4>) ~ Acos(0). 
For t/T > 1 gaps are opened again (suggesting that 
Itoti(t>) diminishes, see below). 

Formally, the supercurrent, for A <^ Tk, is still given 
by the expression 

/(0) =Tsin20[l-rsin2(/)]"^^^ (7) 

but now the transmission T is T = 4pf/[((eo — i)^ + 
f 2)((eo + if + f 2)]. We can interpret that for small t/T 
the Cooper pairs hop directly between the two Kondo 
resonances and T <x t^ . The supercurrent presents 
a sinusoidal-like behavior as shown in Fig. El With 
t/T increasing, the physical situation changes drastically 
around t/T = 1, where the Kondo singularities of each 
dot hybridize into a correlated state as a result of the 
coherent superposition of both Kondo states. We find 
T « 1 and consequently the supercurrent-phase relation 
exhibits a typical baUistic-j unction behavior (see Fig.|21l. 
Further increasing t/T makes Itot{4>) smaller, which is 
attributed to the formation of bonding and antibonding 
Kondo resonances. This results in a nonmonotonous be- 
havior of Ic as a function of t/T (shown in Fig. 01 for 
A/T]^ 0.1,0.25,0.5, from top to bottom), with a max- 
imum 'At t — T (coherent superposition of both Kondo 
resonances). For the lowest gap (in the short junction 
regime) the maximum critical current reaches the uni- 
versal value of 2eA//i as expected |2(Tj |. 

The physics of the tunability of Ic as a function of 
t/T is similar to the transistor-like control of supercur- 
rents in a carbon nanotube quantum dot connected to 
superconducting reservoirs recently reported by Jarillo et 
al, [i^l- These experiments demonstrate that the super- 
current flowing through the QD can be varied by means 
of a gate voltage which tunes on- and off-resonance suc- 
cessive discrete levels of the QD with respect to e^? of the 
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FIG. 4: (Color online). Critical current as a function of t/T 
for A/T^ = 0.1,0.25,0.5 (top to bottom). Upper left in- 
set: comparison with the Josephson critical current given by 
Eq. (dashed line) for A/T^ = 0.1. Lower right inset: 
critical phase at which the maximum current occurs. 



reservoirs. In our case, the Kondo resonances play the 
role of the discrete levels in the experiments of Ref. [2^ 
whereas t/T is the extra knob that tunes the position 
of the Kondo resonances and thus modulates the critical 
current. Interestingly, in our case, i) the supercurrent 
is mediated by a coherent many-body state (the Kondo 
resonance) instead of a single particle one and ii) the 
maximum supercurrent dA, t — T corresponds to coher- 
ent transport of Cooper pairs through the whole device 
whereas an increase oi t/T > 1 splits the Kondo reso- 
nance into two (bonding and antibonding) resulting in 
a splitting of the Cooper pair into two electrons (one on 
each resonance) and thus a reduction of the supercurrent. 
This suggest the use of the Kondo effect as an alternative 
to previous proposals using DQDs for generating and 
manipulating entangled pairs in a controlled way. 

For the experimental realization of the superconduct- 
ing Kondo DQD, we propose carbon nanotubes since (i) 
they show Kondo physics [l^, (ii) it is possible to fabri- 
cate tunable double quantum dots ^2^, and (iii) they are 
ideal systems to attach new material as electrodes . 

Conclusions. — In closing we have studied Cooper pair 
transport through an artificial Kondo molecule. We find 
remarkable differences in the phase-current relation when 
the DQD is built in series or in parallel. For a DQD 
in parallel, the supercurrent always decreases with t/T 
whereas for a serial Kondo molecule the current behaves 
nonmonotonously. This fact allows an extra control of 
the critical current, and thus of Cooper pairs, through 
Kondo molecules by simply tuning the interdot tunneling 
coupling. 
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